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SUPPLEMENTARY NOTE ON MODIFIED -IMPACT -THEORY CALCULATIONS 
FOR BODIES OF REVOLUTION HAVING MINIMUM DRAG AT 
HYPERSONIC SPEEDS 
By Meyer M. Resnikoff 

SUMMARY 


Following the methods of NACA RM A51K27, 1952, the modified impact 
theory developed therein is employed to obtain improved expressions for 
calculating the shapes of bodies of revolution having minimum pressure 
foredrag in hypersonic flight, corresponding to cases where the hyper- 
sonic similarity parameter (the ratio of the free -stream Mach number to 
the slenderness ratio) is appreciably greater than 1. The investigation 
is carried out for various combinations of the conditions of given body 
length, base diameter, surface area, and volume. A minimum drag body of 
given base diameter and surface area is calculated and compared with the 
cone, the corresponding body obtained with impact theory. It is found- 
that consideration of centrifugal forces in the disturbed flow (with the 
modified impact theory) yields a shape of increased bluntness in the 
region of the nose and increased curvature in the region downstream of 
the nose, which result is in substantial agreement with that already 
obtained in NACA RM A51K27 for the minimum drag body of given fineness 
ratio. The calculated pressure drags at hypersonic speeds of the bodies 
obtained with modified impact theory were only slightly less (of the 
order of a few percent) than those of the corresponding bodies obtained 
with the impact theory. 


INTRODUCTION 


Approximate shapes of nonlifting bodies of revolution having ' 

minimum pressure foredrag at high supersonic airspeeds were calculated 
in reference 1 on the basis 'of Newtonian impact theory. The investiga- 
tion was carried out for various combinations of the conditions of given 
body length, base diameter, surface area, and volume. Comparison between 
theory and experiment indicated that bodies so calculated do indeed have 
relatively low drag; however, it was also suggested in reference 1 that 
centrifugal forces in the disturbed flow about such bodies may 
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significantly influence their shape, particularly if the value of the 
hypersonic similarity parameter is appreciably greater than 1 (i.e., at 
hypersonic flight speeds). The impact theory was therefore modified to 
account for effects of these forces in hypersonic flight and the corre- 
sponding minimum drag body of given length and base diameter was calcu- 
lated. This body had more bluntness in the region of the nose, and 
more curvature in the region downstream of the nose than the shape 
obtained with the impact theory. 

Since minimum drag bodies obtained with the modified impact theory 
would appear to be especially suited for hypersonic flight, it has been 
undertaken in the present report to extend the calculations of refer- 
ence 1, using procedures paralleling those presented therein. In 
particular, the modified theory is employed to develop expressions for 
calculating minimum drag bodies of given base diameter and body surface 
area or volume, and given length and volume or surface area. 


SYMBOLS 


ID 


Ctj drag coefficient 2 

Lqo«(*/2) 2 J 

D pressure foredrag 
d maximum body diameter 
f integrand function 


drag parameter ( — — — ) 

V2nq o y 

l body length 

/P-P D \ 

P pressure coefficient ^ J 

p static pressure 

q. dynamic pressure 
S body surface area 


V 


body volume 
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x,y coordinates of point on meridian curve of body (origin of 

coordinate system coincides vith nose of body, and x axis 
coincides with axis of symmetry) 

6 angle (in meridian plane) between free -stream direction and 
tangent to body surface 

X Lagrange multiplier 

Subscripts 

0 free -stream conditions 

1 values at nose point of minimizing curve 

2 values at base point of minimizing curve 

+ right-hand limiting value of quantity at corner on minimizing 
curve 

- left-hand limiting value of quantity at corner on minimizing 
curve 


PROCEDURE FOR CALCULATING MINIMUM DRAG BODIES 
Modified Drag Theory 

It was shown in reference 1 that the pressure coefficient obtained 
by considering only the impact forces may be modified as follows to 
include centrifugal -force effects in the flow over the surface of a 
body: 

t - 2 sin* 8 + | (l - i) £ si* 2 8 (1) 


where the first term on the right represents the contribution due to 
impact forces and the second term represents the effect of curvature 
of the body in the stream direction. A drag parameter may be defined 
by the relation 


Id 



P yy* dx 
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or, considering equation (l), this expression may he written in the 
form (see reference l) 


I D = y a * - f 1 fi + ^-)-^d* 
2 j 0 \ 2y 2 / 1 +y’ 2 


( 2 ) 


Variational methods will he used to minimize equation (2), subject to 
the following conditions: 

a. Given hase diameter and given surface area 
h. Given hase diameter and given volume 

c. Given length and given volume 

d. Given length and given surface area 


Minimizing Procedure 

Relations for volume and surface area can he written 




y 2 dx = const. 


(3) 


and 


A- = + J l yj 1+y' 2 dx = const. 


( 4 ) 


respectively. The conditions for a fixed volume or surface area may he 
included in the expression to he minimized hy adding a multiple of V/fl 
or S/2n to Id, thus forming a new function, say Jp, to he minimized. 
Hence , 


pl 

A + Jy_ N 

y ’ + Xvl 


V 1 

j - + *y 

' 1+y’ 2 j 


( 5 ) 


for a "body with a given volume and 


j d = id + >^= y 2 2 + * 


K 

2 


+ r l + ^ + xyi+y* 2 i y ax 

4> L \ 2y a / 1+y' 2 J (6) 


Ci 
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for a tody with a given surface area, where X is a constant which will 
he determined later. 

In order to find the forms of the function y = y(x) which will 
minimize expressions (5) and (6), it is necessary to find solutions to 
f the Euler differential equation (see reference 2) 


_d_ 

dx 




‘ fy = 0 


(7) 


where f = f (y, y* ) represents the integrand function in equation (5) 

or (6) , and f y t and f y denote, respectively, f (y, y* ) 

B • °y’ 

and — f (y, y’ ) . It can he readily verified that a first integral 
dy 

to equation (7) is given hy 


y* fyt - f = const. 


(8) 


'It follows that a first integral to the Euler equation for the given 
, volume condition is, then, 



2y ,a 

(i+y 2 ) 2 


Xj 


y = c 


(9) 


and for the given surface-area condition. 



3y ^ 2y» 3 

2y 2 J (i +y » 2 ) 2 



( 10 ) 


Solutions to equations (9) and (10), satisfying terminal conditions, 
represent meridian sections of minimum drag bodies (excluding any finite 
section of infinite slope at the nose) for the given conditions, and will 
he called minimizing curves. 

Since the ordinate of the minimizing curve at the nose of the body 
is not specified, a terminal condition must be satisfied (reference 3). ' 
For the given volume condition this is 
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and for the given surface -area condition 


0 .[ fy. ±.?L) -y, r -ML + ( 1 + Jl) ZfA X 

y ^ 2 4o 1 L v - ' -=' 2 


2y 2 / (l+y ,2 )‘ 


y=yi 


( 12 ) 

Similarly, a terminal condition at the base must be satisfied, namely. 




(13) 


when the base diameter is not given, and 


(y* V - f ) x =x„ = 0 


(HO 


when the body length is not given. 

Condition (13) is, for a given volume, 


y 2 


5 y ,g -i 

L 2 (l+y’ 2 ) 2 


+ 2 


= 0 


(15) 


y=y c 


and for a given surface area, 

= o (16) 

y=y 2 

Equation (l4) simply requires that the constant of integration be zero 
in equation (8), the first integral to Euler's equation. 


iy’ 


+ 2 z 


•*-i 


-/l+y' 2 2 d+y ,2 ) ; 


+ 2 


In addition to the above, two further necessary conditions of the 
calculus of variations are needed in order to show that each of the . 
various combinations of given conditions determines a unique min imi zi n g 
curve, namely, the Legendre condition (see reference 2), which requires 


f y ' y ' ^ 0 (17) 

everywhere on the minimizing curve, and the corner condition (see refer- 
ence 2) 
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(18) 


which must be satisfied if the minimizing curve is to contain a corner. 
It can be shown that equation (l8) cannot be satisfied by the integrand 
function of expression (5) or (6). It follows, then, that there can be 
no corners on any of the minimizing curves to be considered, except at 
end points. 


CALCULATION OF MINIMUM DRAG BODIES 
Given Base Diameter and Surface Area 


The first integral to Euler's equation and the terminal condition 
at the base (equations (10) and (l4), respectively) permit the minimiz- 
ing curve to be represented parametrically as 


y = 


x = 


y 2 

3 


d+y^) 


2x3/2 


y dy 

r 


y* * 

xy g 

= 8 




_ 3/2 
2(l+y') 

L y’ 4 


J l+y' 


, c 

y' 


- It\ 




l+y'* 

y' 


> (19) 

! 


The condition at the nose (equation (12)) requires j x = Q or y 1 ' = oo. 
For all values of base diameter (and surface area) the former require- 
ment gives the drag parameter (see equation (2)) a smaller value. Thus, 
it is necessary that y a = 0 and equation (19) yields 


gyx* 3 5y a ' 3 

(i+ yi ’ 2 ) 3/2 (1 +v 2 ) 3/2 


(20) ' 


Using equations (19) and (20), the surface-area condition (4) yields 

X = 121.6 (y 2 2 / S) 3 

From equation (20) it is seen that the range of X is 0 to 2. The 
corresponding ranges for the length l and surface area S are » to 
0.32 y 2 and » to 3*93 y 2 2 > respectively. 
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The minimi zing curve for a body of diameter 2 and surface area 31*57 
is shown in figure 1. The meridian section of a cone (the corresponding 
body obtained with the impact theory) of the same base diameter and 
surface area is also shown in the figure. It is seen that the considera- 
tion of centrifugal forces results in a body of greater bluntness in the 
region of the nose and greater curvature in the meridian section down- 
stream of the nose. Calculation of the drag parameters for these bodies, 
u^ing equation (2), indicates that the cone will indeed have the higher 
drag at hypersonic speeds, although not by more than a percent or two. 
This result and the results of corresponding calculations for the given 
fineness -ratio bodies treated in reference 1 indicate that consideration 
of centrifugal forces principally influences the shape and not the drag 
of minimum drag bodies. 


Given Base Diameter and Volume 


The first integral to Euler’s equation and the terminal condition at 
the base (equations (9) and (1*0, respectively) permit the minimizing 
curve to be represented parametrically by the relations 


where 


y = 


x = 


g y 2 

Xy 2 (l+y» g r 

.3 

y’ 

f r 

U^ r J 


- 3 


x y , = 5 y ' 3 /(i + y 2 ,2 )‘ 


(21) 


The nose condition (equation (ll)) requires either that y x ’ ■ 1, 
or y x ' = ®, or y 1 = 0. If y 1 ' = 1, the resulting if d ratio would 
be limited to values less than one -half, since dy/dy* is positive. 

If y x »’ = oo , the Legendre condition would be violated. Thus it is con- 
cluded that, at least for if d ratios greater than one -half, y x = 0. 

It follows in this case from equation (21) that y x ’ = 0. The Legendre 
condition, namely, 
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limits the range of y 2 * to 0 < y 2 ' < The corresponding range 

of X is from 0 to 3 .■/3/16. 


Given Length and Volume 


Using the first integral to Euler’s equation (equation (9)), the 
following parametric^representation is obtained for the minimizing 
curve: 


y = 


x = 


-y |3 + J f' G + 


C (i+y ,2 r 




3y* 


y 2 


-x (l+y* ) 


» 1 2 \ 2 


3i 


dy 

y* 


(22) 


The terminal conditions at the nose and base are given by equations (ll) 
and (15), which yield the values ’ =1 and y 2 * = 0.27^-, respectively, 

(it can be shown that y 1 cannot be 0) . 

The range of X is -» to 0, corresponding to the range 0 to » 
for l. For X = 0 this becomes the given length and diameter case 
(see reference l) . 


The numerical integration of equations (22) is accomplished by 
first evaluating equation (9) at y = y 2 and y 2 ' = 0 . 27 k and solving 
for c/y 2 in terms of y 2 X. Letting cp(y* , y 2 X) represent the resulting 
function of y* and' y £ X, equations (22) now give 

y = y 2 <p(y f , y 2 M 1 


O .274 


„ [ 7 * £ = y 2 r d< P ( y'.,ya^ L , A (y 2 x) 
4i r y ’ 


1 If y were zero, then it would follow that c = 0 in equation (9) • 
Equations (9) and (15) would then prescribe the value of X. Thus the 

problem would be overdetermined, that is, the above conditions could 
not be satisfied for a general body length and volume. 
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and the volume (equation (3)) is given by 



The values of the functions A and V are obtained by numerical integra- 
tion for various (estimated) values of y to enable interpolation for 
that value of y 2 \, which makes 

r_ = VA 

A 3 x 3 


The set (y 2 X, A, r) so determined satisfies the given volume and length 
requirements and yields the base ordinate value 

y 2 = i/a 

Expressions (22) now determine y parametrically as a function of y» 
and, by numerical integration, x as a function of y* . 


Given Length and Surface Area 


The first integral to the Euler equation (equation (10)) enables the 
minimizing curve to be represented parametrically by 


y = 


_ [X(l+y tg ) 3/2 -2y |3 ] y^±J [X(l+y ,a ) 3 ^ a -2y t3 ] a + 12 c y. 3 (i +y ,3) ; 


6y 


» 3 


. P J *L 

J Q y* 


( 23 ) 


The terminal conditions at the nose and base are given by equations (12) 
and (16), respectively. As in the previously discussed case, the problem 
is overdetermined if y x = 0. Hence, the body in this case must have a 
i mite blunt nose. 
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The range of X is given by 

' -0.64 < X < 2 ^1+ 

As X approaches -0.64, the body length approaches zero. 

The procedure used to integrate equations (23) is similar to that used 
to integrate equations (22) above. 


CONCLUDING REMARKS 


The modified impact theory of NACA RM A51K27, 1951, was employed to 
determine improved expressions for calculating the shapes of bodies of 
revolution having minimum pressure foredrag at hypersonic speeds. Vari- 
ous combinations of the conditions of given body length, base diameter, 
surface area, and volume were treated. A minimum drag body was calcu- 
lated for the case of given base diameter and surface area, and was 
found to be blunter in the region of the nose and to have more curvature 
in the region downstream of the nose than the corresponding shape, a 
cone, obtained with impact theory. Centrifugal force effects considered 
by the modified impact theory were found to influence principally the 
shape and not the drag of minimum drag bodies. 


Ames Aeronautical. Laboratory 
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Moffett Field, Calif. 
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